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Abstract. We consider cosmological scenarios originating from a single imperfect fluid with 
bulk viscosity and apply Eckart's and both the full and the truncated Miiller-Israel-Stewart's 
theories as descriptions of the non-equilibrium processes. Our principal objective is to inves- 
tigate if the dynamical properties of Dark Matter and Dark Energy can be described by a 
single viscous fluid and how such description changes when a causal theory (Miiller-Israel- 
Stewart's, both in its full and truncated forms) is taken into account instead of Eckart's 
non-causal theory. To this purpose, we find numerical solutions for the gravitational poten- 
tial and compare its behaviour with the corresponding ACDM case. Eckart's and the full 
causal theory seem to be disfavoured, whereas the truncated theory leads to results similar 
to those of the ACDM model for a bulk viscous speed in the interval 10~^^ ^ < 10~^. 
Tentatively relating such value to a square propagation velocity of the order of T/m of per- 
turbations in a non-relativistic gas of particles with mass m at the epoch of matter-radiation 
equality, this may be compatible with a mass range 0.1 GeV < m ^ 100 GeV. 

Keywords: Dark Matter, Dark Energy, Gravitational potential, bulk viscosity, Eckart's 
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1 Introduction 

On large scales our present universe seems to be dynamically dominated by two so far un- 
known components called Dark Energy (DE) and Dark Matter (DM), at least if the validity 
of Einstein's General Relativity is taken for granted. To uncover the physical nature of these 
components is a major activity within the cosmologists' community. Lacking a fundamental 
theoretical understanding, the attempts to interpret the wealth of observational data has 
led to a host of phenomenological models with varying degrees of physical motivation. The 
most favoured one is the ACDM model, which is also the simplest. Even though this model 
has good chances to be the "correct" one, it continues to make sense to study alternatives 
that may serve to complement our understanding of the cosmological dynamics, let alone the 
cosmological constant problem [1]. 

Given the large variety of phenomenological approaches, it seems desirable to look for 
unifying aspects of at least some of them. This may also imply a modified view on the status 
of the ACDM model. For example, the latter can be regarded as a particular case of the 
generalised Chaplygin gas (GCG) models [2-5]. In turn, it was shown that the GCG can be 
interpreted on the basis of interacting holographic DE models [6, 7], where the interaction 
parameter in the dark sector determines the GCG parameters. 

A further line of research have been dissipative DE models, in which the negative pres- 
sure, responsible for the current acceleration, is an effective bulk viscous pressure, i.e. a 
non-equilibrium phenomenon. This parallels related approaches for an inflationary phase in 
the early universe [8-12]. It was argued in [13, 14], that a viscous pressure can play the role 
of an agent that drives the present acceleration of the Universe. Note that the possibility 
of a viscosity dominated late epoch of the Universe with accelerated expansion was already 
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mentioned by Padmanabhan and Chitre in 1987 [15]. For more recent investigations see, e.g. 
[16-21] and references therein. It has also been pointed out in [16-18] that there exists a 
correspondence between dissipative and GCG models such that for appropriate parameters 
combinations both models result in an identical dynamics. 

It is expedient to emphasise, however, that these very different manners of understanding 
the apparently same dynamics, say of the ACDM model, are restricted to a homogeneous and 
isotropic background. At the perturbative level, the different mentioned approaches imply a 
generally different dynamics. The unifying view on the background level is accompanied by a 
difference in the perturbation dynamics, e.g., as a consequence of fluctuations in the relevant 
interaction rates [6, 7]. This circumstance makes the different approaches potentially testable 
alternatives to the "standard" ACDM model while sharing the same background dynamics 
of the latter. 

Traditionally, for the description of non-equilibrium thermodynamical processes the the- 
ories of Eckart [22] and Landau and Lifshitz [23] were used. Due to the work of Miiller [24], 
Israel [25], Israel and Stewart [26], Pavon, Jou and Casas- Vazquez [27], Hiscock and Lind- 
blom [28] it became clear, however, that Eckart-type theories suffer from serious drawbacks 
concerning causality and stability. These difficulties could be traced back to their restriction 
to first-order deviations from equilibrium. If one includes second-order deviations as well, the 
corresponding problems disappear. Cosmological implications of second-order theories were 
first considered by Belinskii, Nikomarov and Khalatnikov [29], followed by Pavon, Jou and 
Casas- Vazquez [27] and other authors [9-11]. Reheating and cosmological particle production 
on the basis of causal thermodynamics have been discussed in [12, 30]. 

One should be aware that discussing the issue of bulk- viscous driven accelerated expan- 
sion at all, implies in any case an extrapolation of non-equilibrium thermodynamical theories 
beyond the range for which their applicability was strictly justified. Bulk viscous infiation, 
if it exists, is a far-from-equilibrium phenomenon, while even the full causal second-order 
Miiller-Israel-Stewart's (MIS) theory holds for small deviations from equilibrium. There- 
fore, all theoretical conclusions are necessarily tentative but, hopefully, they nevertheless will 
provide an indication of the correct behaviour far from equilibrium. 

The impact of bulk viscosity on the background expansion of the universe is widely 
investigated in the literature, see e.g. [9, 18, 29, 31-42]. However, the perturbative analysis 
of the viscous cosmological models is not as widely addressed, despite its crucial importance. 
Among those papers which do perform a perturbative analysis, we cite [16, 17, 20, 43, 44]. 
A common characteristic of these papers is the analysis of a fiuid whose equilibrium pressure 
is negligible with respect to the bulk viscous one. Moreover, another common point is that 
the dissipation process is described within Eckart's theory. Interestingly, there exist choices 
of the bulk viscosity parameters which give predictions for the matter power spectrum in 
agreement with observation [43, 44]. On the other hand, the results found in [20] seem to 
doom the possibility of describing DM and DE together as a bulk viscosity effect because a 
very strong integrated Sachs- Wolfe effect would be produced. 

The present paper relies on a cosmological dynamics, realized by a single imperfect 
fluid with bulk viscosity and vanishing equilibrium pressure. However, we investigate how 
the picture changes when taking into account the causal Miiller-Israel-Stewart theory, both 
in its full and truncated forms, instead of Eckart's theory. Our focus is on the behaviour of 
the gravitational potential, which we compare with its ACDM counterpart. Assuming that 
the results of the standard ACDM model grosso modo reproduce the observations, we will 
discard any viscous model that differs substantially from this standard. 
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For both Eckart's and the full MIS theories, our results seem to disfavour the possibility 
of unifying DM and DE via bulk viscosity. Only the truncated case seems to remain a valid 
option for bulk viscous speeds of sound < 10~^. It may be interesting to mention that 
extremely small values for the speed of sound are also known from certain scalar-field DE 
models [45]. 

The paper is structured as follows: in Sec. 2 we present the fundamental equations of 
bulk viscous cosmology. In Sec. 3 we set up the formalism necessary to investigate pertur- 
bations in the bulk viscous fluid. In Sec. 4 we adopt and motivate an assumption for the 
viscosity pressure and in Sec. 5 we show our results, embodied in the evolution of the gravi- 
tational potential for the viscous models compared to the ACDM one. In Sec. 6 we draw our 
conclusions. We choose the —2 signature for the metric and c = 1 units. 



2 Dynamics of the viscous fluid 

A flat, homogeneous and isotropic universe is described by the Robertson- Walker (RW) line 
element 

ds'^ = dt^ - a{tf5ijdx'dx^ = a{r]f {dif - dijdx'dx^) , (2.1) 

where t is the cosmic time and rj the conformal one. Consider a relativistic fluid with bulk 
viscosity. We refer the reader to [9-11] for extensive reviews of the topic. Write the total 
pressure as P = p + 11, i.e. consider a splitting into the equilibrium part p plus the bulk 
viscous pressure contribution 11. In presence of the latter, the stress-energy tensor has the 
following form: 

T^j,y = {p + P) u^Uy - pg^y , (2.2) 

where is the four- velocity and gfj,^ is given by Eq. (2.1). Hence, the bulk viscous pressure 
part simply adds up to its equilibrium counterpart. 

Neglecting the baryon and the radiation components,^ Friedmann's equation and the 
energy conservation equation read 

i/2 = ^ = ^p, p = -2,H{p + P), (2.3) 



a 

where H = a/a \s the Hubble parameter and % = a' /a. Throughout the paper, the dot 
denotes derivation with respect to (wrt) the cosmic time whereas the prime denotes the 
derivation wrt the conformal time. Combining the two equations in (2.3), one obtains the 
acceleration equation 

which can be solved for P/ p, resulting in 

^ = ^(2^-1), ^'-=-^-§2^ (2-5) 

where q is the deceleration parameter. 

^The baryon component is always subdominant whereas radiation is negligible well after decoupling, which 
is the timespan we are interested in. 
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The evolution of 11 is governed by a transport equation whose form depends on the 
theory adopted for the dissipative processes. Eckart's theory [22] is based on the algebraic 
relation 

u = -ec, (2.6) 

where C is the bulk viscosity coefficient and 9 = V^u^ is the expansion scalar, i.e. the 
covariant divergence of the velocity field. For the RW metric one has 6 = 3H = 3T-L/a. On 
the other hand, in the MIS theory [24-26] the evolution of IT is described by the following 
transport equation: 



rn* + n 



1 



-rn 



(C/r)' 

(CM 



(2.7) 



where r is a relaxation time and T is the temperature. The bullet • in Eq. (2.7) denotes 



n* 



(2. 



i.e. derivation along the fluid wordline. In the RW case it reduces to the derivation wrt 
the cosmic time, but we prefer to keep a separate notation in order to avoid confusion when 
considering perturbations of (2.7). Equation (2.7) has often been employed in a truncated 
form, i.e. 

rn' + n = -0c, (2.9) 

which we are going to investigate in detail, comparing it with the full theory and Eckart's 
one. 

Note that, assuming the bulk viscous pressure 11 to be subject to causal thermodynamics 
does not only introduce the relaxation time as an additional parameter but also, associated 
with the relaxation time, a propagation velocity for viscous perturbations which is different 
from the adiabatic sound speed. 

The temperature term in Eq. (2.7) can be calculated from Gibbs' integrability condition 

[11]: 

dT , ^dT ^dp 



^ dn~^ ~^ dp '^dp' 



(2.10) 



Assuming an equation of state p = p{p, n), integration of (2.10) provides us with a functional 
form for T = T{p, n). Then, energy and particle number conservations allow us to determine 



dp n dT 
dp ^ T~dp 



(2.11) 



In this paper, we assume the temperature and the equilibrium pressure to be barotropic, i.e. 
described by T = T(p) and p = p{p). Defining dp/ dp = dp /dp := Cg as the adiabatic speed 
of sound, Eq. (2.10) takes the form 



and Eq. (2.11) becomes 



1 dT 
T dp 



T' 

^ = -Ocl 

rj-i S 



p + p 



1 + 



n 



p + p 



(2.12) 



(2.13) 
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(With our assumptions T = T{p) and p = p{p) we disregard a possible microscopic description 
on the basis of kinetic theory). Under this condition, we can cast the MIS transport equation 
in the following form: 



rn' + n 



-9C 



1 



-rn 



{CM 



1 



n 



p + p 



(2.14) 



The coefficients r and Q are, in general, functions of the time (or of the energy density). They 
are not completely arbitrary but related to c^, the square of the velocity at which a viscous 
perturbation propagates [28] (see also [11]): 



c 



{p + p)t 



(2.15) 



Remarkably, this propagation speed depends also on the equilibrium pressure contribution. 
Given that has to be added to Cg, for causality reasons we have 



(2.16) 



which provides a constraint on a possible functional form of r and C, once the adiabatic speed 
of sound is given. 



3 Perturbative dynamics of the viscous fluid 

In this section we provide at first the general perturbation equations which then are specified 
to Eckart's and MIS theories. We adopt the Bardeen gauge-invariant formalism [46-48]. See 
also [49] for an extensive review of Eckart's case. 

Consider the following perturbation of the RW metric: 

ds^ = a{r]f [(1 + 2$)dr/2 - (1 - 2^)d^^] . (3.1) 

The background velocity field has components uq = a and Uj = 0. The perturbed zero- 
component is 6uo = a^, 6u^ = — $/a. A general perturbation of the stress-energy tensor 
(2.2) has the form 

6T% = Sp , (3.2) 
= ip + p + U)v\ (3.3) 
ST'j = -6'j {6p + <5n) , (3.4) 

where we have defined := a5u^. At first order, Einstein's equations take the form 

A$ - 3n ($' + n^) = A7rGa^5p , (3.5) 

A{<^' + n<^) + {n^ -n')e = , (3.6) 

+ 3n^' + (^2 + 2n') $ = 47rGa2 (Sp + 5U) , (3.7) 

where we have used the divergence of the (0 — i) equation and defined © = div''. The 
expression for the perturbed bulk viscous contribution 511 depends on the choice of the 
theory. 
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3.1 Eckart's theory 

From Eq. (2.6), a perturbation of IT may be formally written as 

6u = -6ec - esc ■ 

With e = Vf,u^' we find 



(3.8) 



56 = d^6u'' + ri'p^duP + Sri'p^uP , (3.9) 

where Fp^ is the Levi-Civita connection of the RW metric and STp^ is its perturbation. Start- 
ing from Eq. (2.1) and Eq. (3.1) and working out the above expression, it is straightforward 
to obtain 

a69 = dy - 3 ($' + n<^) = e - 3 ($' + n^) . (3.10) 

Therefore, 

3^ _ C 



du = -—6c - - [G - 3 (^>' + n^)] , 

a a 



(3.11) 



is the perturbed bulk viscous pressure contribution we are looking for. 
3.2 Miiller-Israel-Stew^art's theory 

In the MIS theory, the bulk viscous pressure contribution evolves according to Eq. (2.14). In 
the following we shall always replace C by Cb. Note that this is not possible for Eckart's case 
since it has r = 0, so that diverges, cf. Eq. (2.15). In order to simplify the notation, let 
/ = /(/?) := C/t. This choice also leaves open the possibility of a variable c^, although we 
do not consider this case in the present paper. 

Equation (2.14) can then be cast in the compact form 



r 



where we have defined 



gip) :=2 + c2(p) + 



P + P{p) dclip) 



fip) + ^9{P) 



Hp) 



Hp) 



ctip) dp 

The derivative term in Eq. (3.12) is perturbed as follows: 

6 (W) = 6 (u'^d.U) = --U' + -6U' 

a a 

The first-order equation that results from Eq. (3.12) is 



1 dclip) ^ 1 + 2cUp) 
cUp) dp p + p{p) 



(3.12) 



(3.13) 



(3.14) 



1 



n 



-6U' + ,5 - = -n' 



fip) + ^9ip) + ^Hp) 



dfip) ^Ildgipl^U^ dHp) 



9511 



9{P) 



dp 2 dp 2 dp 
Substituting Il'/a from the zeroth order of (3.12), one obtains 
1 



+ Uhip) 



(3.15) 



-5U' + 6 



a 



6^ 



fip) + ^9ip) + Y^ip) 



esp 



dfjp) ^ Udgjp) ^ nUhjp) 



dp 



2 dp 



dp 



-56 
6511 



fip) + ^9ip) + ^Hp) 



9ip) 



+ Uh{p) 



(3.16) 
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This equation determines the perturbation 511 that enters the right-hand side of (3.7). 
The truncated form of (3.16) can be obtained by putting g = h = 0: 

-SW + 5('^)=-^- emp) - mip) - OSp'Jj^ . (3.17) 
a \t J T dp 

In the fohowing section we analyse special forms of C and r which provide a background 
expansion similar to the one produced by the generahsed Chaplygin gas [2, 3]. 



4 Ansatz for the pressure 

Up to now, our formalism for a one-component bulk viscous fluid with barotropic equation 
of state has been quite general. Now and hereafter we also assume the effective pressure to 
be totally dissipative, i.e. P = U. This implies p = and = 0. Moreover, we demand 
to be a constant, so that in MIS theories we have C = '^hP'^- 

In the causal theory the viscous pressure becomes a dynamical degree of freedom, subject 
to a first order differential equations, both in the background and on the level of perturba- 
tions. However, in the context of this paper these differential equations are not used to 
determine the viscous pressure and its perturbation. Instead, we shall assume the viscous 
pressure to be known. Then, the causal evolution equation provides us with a relation be- 
tween the relaxation time and the Hubble time. This strategy is motivated as follows. Our 
present Universe is well described by the ACDM model, both at the RW level and for linear 
perturbations about that background. As far as the former is concerned, the ACDM model 
corresponds to a constant pressure in a one-component model of the cosmic substratum. In 
order to admit also deviations from the ACDM model, we start with an ansatz 

(4.1) 

P \PoJ 

where p is a, constant and pQ is the present time (i.e. for a = a^) energy density of the 
viscous fluid. For n = 2 we recover the constant pressure case. Equation (4.1) represents the 
equation of state of a generalised Chaplygin gas [2, 3]. Consistently, by solving the energy 
conservation equation one can find, taking oq = 1, 



P = Po 



1 2/n 

p + {l-p) a-3"/2 . (4.2) 



In the limit a <^ 1 the solution (4.2) reproduces a matter dominated universe with p oc a~^, 
while in the opposite limit the energy density is similar to that of a cosmological constant. 
Moreover, for n = 2 we recover the ACDM model while for n = A the expression (4.2) 
describes the energy density of the "true" Chaplygin gas. 

Plugging (4.1) and (4.2) into (2.4) and (2.5), the deceleration parameter can be written 

as 



1 

'^=2 



1- - ^ 



p + {l- p) a~'^'^l'^ 



(4.3) 



Denoting the value of a at g = by Oq, we have 
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In terms of either Oq or qo the pressure to energy density ratio takes the form 
n 1 1 - 2go 



Oq \ 3n./2 



and hkewise the Hubble rate becomes 



1 - 2^0 + 2 (1 + qo) a-3"/2 ' 



(4.5) 



H 



l + 2(^)^"/2 
l + 2ai^"/2 



l/n 



1/n 



1 - 2^0 + 2 (1 + qo) a 



-3n/2 



l/n 



(4.6) 



The parameters ^ = C(p) ^ind r = r(p) describe intrinsic properties of the fluid and 
their perturbations have the following forms: 



dp 



dp 



(4.7) 



Equation (4.1) then provides us with the structure for Cip)- fhe causal theories it will 
additionally allow us to calculate the relaxation time t{p). 

4.1 Eckart's theory 

From (4.1) and H = —9(^ we can solve for and find 

^l-n/2 



C = PPq 



lip 



n/2 

PPo 

p ^ 



3H V24^ 



(4.8) 



where we have used Friedmann's equation (2.3) to express H as function of p. Taking into 
account (4.7) we find then 

f = ii^^ . (4.9) 

C 2 p 

which we are going to use in the perturbative analysis of the next section. 
4.2 Full Miiller-Israel-Stewart's theory 

In this section the • derivative reduces to the time derivative and will be denoted by a dot 
again. Assuming to be constant, we find from (3.13) 



g = 2 , and h 



1 

P 



(4.10) 



With n/p given, equation (3.12) provides us with an expression for the parameter tH: 



3tH 



n 

' p 



2 • 



(4.11) 



Using now (4.1) and (2.3) we find for t{p): 



rip) 



2p 



-n/2 



Po 



(4.12) 
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Taking into account (4.7) we obtain 



6t n 6p 



2cl-fi^{n-l) 



~n/2 



1 6p 

2 7 



(4.13) 



To make physical sense, the parameter tH has to be positive. Combining (4.11) and (4.5), 
the condition tH > is equivalent to 

1 i + 2n(^)'"/' 



l + 2(^)^"/^ 



o (1 - 2<7o) 



1 - 2go + 2n (1 + qp) g-^"/^ 
[l-2go + 2(l + go)«"=^"/2]' 



(4.14) 



For arbitrary a, this is only guaranteed for > ^. Remarkably, this limit does not depend 
on n. As we shall see, such a large minimum value for the bulk viscous speed of sound has 
consequences for the perturbation dynamics, which disfavour the theory. 



4.3 Truncated Miiller-Israel-Stewart's theory 

In that case we have 

ri + in = -3Hf , 

r 

instead of (3.12). Solving for the relaxation time we obtain 



3tH 



n 

p 



(4.15) 



(4.16) 



p V 'P. 

Using again (4.1) and (2.3) we find t{p) in the truncated case: 



1 2/.( 




-n/2 


2c^ + M2-n)(^ 


r 


/2 




-n/2" 



(4.17) 



With (4.7) it follows that 
n 6p 



5t 

T 



2c? -/i2(n-2) 



2 P 0.2 



2c^ - /i (n - 2) 



-n/2 



+ /^Mn-2)(^ 



1 6p 

2 7 



(4.18) 



As far as the viscous speed of sound is concerned, the situation is different from that of the 
full theory. Here we have 



3Ht 



1 + 2 (^) 



3n/2 



or, equivalently, 
3HT = il-2qo) 



(^)'"/'] +(2-n) (^)'"/' 
l-2go + 2(l + go)a-3"/2 



(4.19) 



c2 [1 - 2(?o + 2 (1 + go) a-3-/2] V (2 - n) (1 - 2qo) (1 + go) ^-3-/2 



(4.20) 
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In the truncated theory the requirement 3Ht > impUes that 

1 - 2go + 2 (1 + qo) a -^"/^ 

which (for > ) is automatically satisfied for n < 2. In particular, at high redshifts, 
corresponding to a ^ 1, it follows that 

2cl 1 + qo ^ ^ 

This means that either for a ^ 1 and ^ 1 or for <C 1 with ^ a^"/^, one has 

T « //-^ . (4.23) 

A small relaxation time at high redshift is consistent with the limit |n|//3 ^ 1 for a <C 1 [cf. 
(4.5)]. In fact, during this period we have close-to-equilibrium conditions. If we assume our 
model to be valid from the epoch of matter-radiation equality at Zeq ~ 3.3 • 10"^ on, i.e. for 
a > Acq ^ 3 • 10~^, the condition for a model, similar to the ACDM model with n ~ 2, is 
^ 10~^^. As we shall show below, acceptable values for are of the order of < 10"*^, 
which is well compatible with that lower bound. This illustrates that a finite relaxation 
time requires a non- vanishing, albeit small, value of c^. It is obvious from (4.22), that Ht 
increases with the scale factor. 

At the present time a = 1 the expression (4.20) specifies to 

-^^lo = 7 ^ 7 , , r . (4.24) 

'° 9cl + (2 - n) (1 - 2qo) (1 + qo) ^ ^ 

We have 

t|o > . (4.25) 

The relaxation time is of the order of the Hubble time or even exceeds the latter. Finally, in 
the long time limit, i.e., a ^ 1, relation (4.20) becomes 

_ 1 

3 c; 



2 H^' (a > 1) . (4.26) 



It follows, that for any < 1/3, but larger than a very small minimum value ~ 10 the 
relaxation time will remain finite and larger than the Hubble time, i.e., r > H^^ or even 

The situation is different for n > 2. Equation (4.21) is equivalent to the existence of 
a lower limit on c^. More precisely, the maximum of the function on the right hand side of 
(4.21) occurs at 

3n/2 1 

' = 2 • (4-27) 
Therefore, the relevant limit is 

cl>^. (4.28) 

For example, for the Chaplygin gas (n = 4), it amounts to > 1/4. This value is slightly 
smaller than that of the full theory but, as the numerical analysis will show, unless n is 
extremely close to 2, this limit will not admit a favourable scenario. 
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In the next section we shall see that a perturbation theory based on the truncated 
version for n < 2 fares better with respect to the behaviour of the gravitational potential 
than the full theory. Different from the full theory, the truncated theory allows for a finite 
but small value ^ 1. 

It is expedient to recall that both the full and the truncated theories are applied here 
under conditions where a microscopic justification is missing. Both versions were taken here 
on a purely phenomenological basis. Intuitively, the simpler truncated version appears even 
more transparent since it represents the minimally possible generalisation of Eckart's theory. 
In the present context, the truncated theory is not necessarily an approximation to the full 
theory, but it should be considered as a phenomenological description on its own. From this 
point of view one should not be too surprised to obtain "better" results from the truncated 
theory. 

5 Results: the evolution of the gravitational potential 

In this section we numerically solve the equation for the gravitational potential, for the 
Eckart, the MIS and the truncated MIS cases. Independently from the transport theory, the 
background evolution is the same, viz. (4.6), due to our choice (4.1). In particular, for the 
case corresponding to the ACDM model, n = 2, we have 

^2 ^ l + 2(^)^ 



Hi l + 2a3 



(5.1) 



with the ACDM identifications 



1 

For r^Ao = 0.734, we have = 0.566, or a redshift Zq = 0.767 and qq = —0.204. In this 
section we denote the derivative wrt the scale factor by a subscript a. 

5.1 Eckart's case 

Linking 6p to the gravitational potential via the (0 — 0) Einstein equation (3.5), we get from 
(4.9) 

f = ^ [^^ - 3^ + ^^)] • (5-3) 
Using now the (0 — i) Einstein equation (3.6) written as 

= n'-n- ' ^'-'^ 



together with (3.11), (5.3) and the (i — j) Einstein equation (3.7), we obtain 



+ 3n^' + {n'^ + 2n') $ = AnCaC 



(5.5) 



In terms of derivatives wrt the scale factor (subscript a) we get 



(n - 1)A3. - + - mHn - 2) + 



(5.6) 
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Performing a plane wave expansion $ cx exp (zk • x) and using (4.8), we can cast (5.6) in the 
form 



(1 - 2qo) H^k ^ 



(n - 1)^> 



n (g^g + ^) 



(1 - 2qo) 



2Wd 



rin'i—n 



(n - 2) {a<^a + ^) . (5.7) 



For the case that simulates the ACDM model, n = 2, Eq. (5.7) becomes 

The ACDM model is recovered for a vanishing right-hand side. A non-vanishing right-hand 
side is the result of pressure perturbations. It is a characteristic feature of the viscous model, 
that the k"^ term in Eq. (5.7) does not only contain $ but the derivative $a as well. This 
coupling to <I>a (or its energy-density counterpart) turned out to be essential for an adequate 
description of the matter power spectrum [43, 44]. Generally speaking, for the behaviour of 
the gravitational potential such a coupling causes deviations from the ACDM model which 
are especially relevant from the viewpoint of the integrated Sachs- Wolfe (ISW) effect [50], 
where an integration of $ up to /c — )■ oo has to be performed. 

Numerical solutions of Eq. (5.7) are displayed in Fig. 1 and Fig. 2. We plot the ratio 
of the viscous gravitational potential computed from Eq. (5.7) to the ACDM one (which is 
also computed from Eq. (5.7) with k = and n = 2). We have chosen initial conditions 
^(o*) = 1, ^(o*) = 0, where a* = 10"'^ is the recombination scale factor. Moreover, we have 
also chosen = 0.566, as the reference transition redshift of the ACDM model. 

From the plot in Figs. 1 and 2 we can infer the following features: 

1. On the largest scales and for small values of the scale factor, the results of Eckart's 
theory approach those of the ACDM model. On smaller scales and for scale factors 
of the order of one, the differences are substantial or even dramatic. Given, that the 
equations just differ by fc^ terms, this is not surprising. 

2. Varying the background expansion implies the following: for increasing n, the effect 
of the viscosity contribution on intermediate and smaller scales seems to be "delayed" 
and the agreement with the ACDM result enhances for a < Oq. On the other hand, for 
a > the differences increase. 

3. The shape of $/^'a close to a = 1 in Fig. 2 does not depend much on n. An increasing 
n causes only a shift towards smaller values of the ratio $/<I>a- 

Already from this qualitative analysis we may expect the agreement of this model with 
observation to be poor. 
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Figure 1. Ratio $/<I>a in Eckart's theory as function of the scale factor for fc = 5 • 10 10 '^,10 ^ h 
Mpc^^ (soUd, dashed, dot-dashed lines, respectively). From the upper-left panel in clockwise sense, 
n = 1,2,4,8. 
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Figure 2. Ratio ^/^a in Eckart's theory as function of the wavenumber fc for a = 1 fixed. Left 
Panel: n= 1,2 (solid, dashed). Right panel: n = 4,8 (sohd, dashed). 



5.2 Full Miiller-Israel-Stewart's theory 

With our ansatz / = c^p and constant (which imphes g = 2 and h = 1/p), Eq. (3.16) can 
be cast into the following form: 



5U 

P 



1 5Ii 



+ 



1 n(5r 



Har p Har p r 



59 



n 



a \ 2) p 

n2 



Ha + p + 2p2^ ap 2p2 



3 5p 



1 + 5 

n2 



(5.9) 



where the subscript a again denotes derivation wrt the scale factor. The evolution of II/p in 
(5.9) is known by (4.5). Taking into account (4.13) we find the following system of coupled 
equations: 



5Yi\ 3 
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4 + ^1* 
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+ 3(2 - n)^a — 



2p2 
n nX^ 



2p 2p^ 



+ 



+ 



2p 



(2 - n)k^ 



nn (n - i)n2 



2p2 

n 



n 



2p 2p2 



n 



Cb + -: + 



2p^ 



^-n 



an.a — n. 



(5.11) 



Obviously, the simplest case is n = 2, the case which reproduces the ACDM background. It 
eliminates various contributions, among them also a A;2$ term. However, the above system 
of equations is so complicated that only a numerical analysis will uncover its secrets. We 
have now one parameter more than in Eckart's theory, the constant c^. We present numerical 
results of the above system, recalling the limit > 1/2 [see the comment following (4.14)]. 
As for Eckart's case, we choose ^>(a*) = 1, ^{a^,) = as initial conditions on the gravitational 
potential plus 5Ii{a^) = 0, i.e. assuming that perturbations in the viscous pressure of the 
fluid were not important at recombination. Again, we choose Oq = 0.566 as the transition 
redshift. We plot the results for different choices of the background evolution, i.e. for various 
n's, and for different choices of > 1/2. 

We infer the following from Figures 3, 4 and 5: 

1. As for Eckart's case, the ratio <1>/$a decays faster and earlier in time for larger values 
of A:. 

2. For fixed k and increasing n the decay phase of the evolution is more pronounced and 
takes place earlier than in Eckart's case. 



3. 



Apparently, the value of has a very little influence on the evolution of the gravita- 
tional potential. This is probably due to the fact that must be larger than one-half, 
therefore already too large for allowing for a perturbative evolution close to the ACDM 
case. For n = 4, a decreasing even enlarges the difference between <I> and <I>a. 

In conclusion, it seems that taking into account the full causal theory does not amend the 
problems already encountered in Eckart's case. 
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Figure 3. Results for the case n = 1: ratio $/$a in the full MIS theory for f{p) = c^p with 
constant. Top-left panel: ~ 1. Top-right panel: = 0.7. Bottom-left panel: = 0.5. In the latter 
three panels the profiles are functions of a for k = 0.0001,0.001,0.01 h Mpc^^ fixed (solid, dashed, 
dot-dashed, respectively). Bottom-right panel: dependence on k for a ~ 1; here the solid, dashed, 
dot-dashed lines correspond to the cases cl = 1, 0.7, 0.5, respectively. 



5.3 Truncated Miiller-Israel-Stewart's theory 

With our ansatz / = c^p and constant, Eq. (3.17) can be cast into the following form: 
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where a subscript a denotes a derivative wrt the scale factor. Together with (4.18), we find 
the following system of coupled equations: 
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Figure 4. Results for the case n = 2, i.e. a ACDM background: ratio $/<i>A in tiic full MIS 



theory for /(p) = c^p with constant. Top-left panel: 



1. Top-right panel: = 0.7. Bottom-left 



panel: = 0.5. In the latter three panels the profiles are functions of a for fc = 0.0001,0.001,0.01 
h Mpc~^ fixed (solid, dashed, dot-dashed, respectively). Bottom-right panel: dependence on k for 
a = 1; here the solid, dashed, dot-dashed lines correspond to the cases c? = 1,0.7,0.5, respectively. 



Also here, the case n = 
simplest. Introducing 11 



2, corresponding to a ACDM background, turns out to be the 
= —fJ-Po constant and the Friedmann equation, we obtain 
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(5.16) 



where we have separated the terms. In the limit = the right-hand side of (5.15) rapidly 
decays and we recover the ACDM dynamics. However, = is equivalent to / = C = 0, i.e., 
such a case does no longer describe a bulk viscous fluid. This corresponds to the circumstance 
that there also exists a lower limit for [see below (4.23)], derived under conditions where 
the credibility of the theory is highest. Consequently, what is needed, is a very small but 
non vanishing value of c^. 

Let us consider how small has to be in order to have a viable truncated viscous theory. 
We choose <^'(a*) = 1, $(a*) = and 5Yi{a^) = as initial conditions, and Oq = 0.566. We 
plot the results for different choices of the background evolution, i.e. for various n's, and for 
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Figure 5. Results for the case n = 4: ratio $/$a in the full MIS theory for f{p) = c^p with 
constant. Top-left panel: ~ 1. Top-right panel: = 0.7. Bottom-left panel: = 0.5. In the latter 
three panels the profiles are functions of a for k = 0.0001,0.001,0.01 h Mpc^^ fixed (solid, dashed, 
dot-dashed, respectively). Bottom-right panel: dependence on k for a ~ 1; here the solid, dashed, 
dot-dashed lines correspond to the cases = 1, 0.7, 0.5, respectively. 



different choices of c^. 

Figs. 6 and 7 show that: 

1. For n = 4, lowering from unity to the limiting value = 1/4 worsens the situation. 

2. From Fig. 6 we find that for k = 0.01 h Mpc"^, a scale which lays at the border of the 
linear regime, the present value of <I>/<I>a is about 0.6. This difference is compatible with 
the uncertainty due to the cosmic variance [48]. Therefore, we may infer the qualitative 
constraint < 10"^. 

In order to understand a value < 10"^ at least roughly and tentatively, we may 
resort to a kinetic picture. According to a kinetic-theory based close-to-equilibrium analysis 
by Israel and Stewart [26], the propagation speed of viscous pulses is of the order T/m <^ 1 for 
nonrelativistic particles of mass m, i.e., of the order of the adiabatic speed of sound. Assuming 
here the radiation temperature at matter-radiation equality Tcq ~ 10^ A' ~ 0.86 • 10~^GeV 
as an upper limit for the matter temperature, a value ~ 10^^ corresponds to a lower 
mass limit at 0.1 GeV. On the other hand, we must also take into account the lower limit 
cl > 10"" [cf. the discussion after Eq. (4.23)]. The range 10"^^ <■ cl < 10"^ then 
corresponds to a mass interval 0.1 GeV ^ m <^ 100 GeV. Note that this rough estimate is 
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Figure 6. Results for the case n = 2, i.e. a ACDM background: ratio ^/^a in the trun- 
cated MIS theory for f{p) = c^p with constant. Top-left panel: = 1. Top-right panel: 
= 10^^. Bottom-left panel: = 10~*. In the latter three panels the profiles are functions of 
a for A: = 0.0001,0.001,0.01 h Mpc^^ fixed (solid, dashed, dot-dashed, respectively). Bottom-right 
panel: evolution in function of k and for a = 1; here the solid, dashed, dot-dashed lines correspond to 
the cases = 1, 10"'^, 10~®, respectively. 



compatible with the predicted mass of DM candidates such as the WIMPs, and in particular 
the neutralino [51]. 
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Figure 7. Results for the case n = 4: ratio $/$a in tlie truncated MIS theory for f{p) = c^p 
witli constant. Top-left panel: = 1. Top-right panel: = 0.5. Bottom-left panel: = 0.25. In 
the latter three panels the profiles are functions of a for k = 0.0001, 0.001, 0.01 h Mpc~^ fixed (solid, 
dashed, dot-dashed, respectively). Bottom-right panel: evolution in function of k and for a = 1; here 
the solid, dashed, dot-dashed lines correspond to the cases = 1,0.5,0.25, respectively. 



6 Conclusions 

We have investigated cosmological scenarios in which the matter content of the Universe is 
given by a bulk viscous fluid. The bulk viscous pressure was described by Eckart's theory 
as well as by the causal Miiller-Israel-Stewart's theory, both in its full and truncated forms. 
The causal theories introduce a relaxation time as an additional parameter into the cosmo- 
logical dynamics together with a finite viscous speed of sound. Imposing a homogeneous and 
isotropic background dynamics with the ACDM model as a special case, the causal trans- 
port equations for the viscous pressure reduce to relations between the propagation speed of 
viscous pulses and the relaxation time. The requirement of a finite relaxation time implies 
upper and lower limits for the viscous sound speed. We calculated numerically the first-order 
gravitational potential for each of the viscous theories and compared it with the ACDM 
result. 

While on the largest scales and for high redshifts all theories are similar, there are 
substantial and even dramatic differences on smaller scales and low redshifts. Since the 
gravitational potential of the ACDM model is supposed to provide an approximately reliable 
"standard" description of the CMB power spectrum, any theory with strong deviations from 
this standard is very likely also disfavoured by the data. On this basis, Eckart's theory 
seems to be ruled out, despite of the fact that it is able to reproduce the matter power 
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spectrum [43, 44]. This result confirms an analysis of the integrated Sachs- Wolfe effect in 
[20]. Likewise, the potentials of the full Miiller-Israel-Stewart's theory differ too strongly 
from its ACDM counterparts to be a promising alternative. The reason here is a lower limit 
1/2 for the square of the viscous propagation speed, which is much too high for a behaviour 
similar to that of the ACDM model. The only potentially valid option is the truncated causal 
theory for propagation speeds within the range 10~^^ ^ ^ 10~^. In the present context, 
the truncated theory is not an approximation to the full theory but a phenomenological 
description on its own. Tentatively, values of the order of 10~^^ ^ ^ 10"*^ can be related 
to a nonrelativistic particle with a mass in the interval 0.1 GeV < m < 100 GeV which 
corresponds to a typical WIMP mass scale. A more quantitative analysis of the CMB power 
spectrum for the causal viscous model will be the subject of future work. 
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